Crossover in the Slow Decay of Dynamic Correlations in the Lorentz Model 
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The long-time behavior of transport coefficients in a model for spatially heterogeneous media in two and three 
dimensions is investigated by Molecular Dynamics simulations. The behavior of the velocity auto-correlation 
function is rationalized in terms of a competition of the critical relaxation due to the underlying percolation 
transition and the hydrodynamic power-law anomalies. In two dimensions and in the absence of a diffusive 
mode, another power law anomaly due to trapping is found with an exponent —3 instead of —2. Further, the 
logarithmic divergence of the Burnett coefficient is corroborated in the dilute limit; at finite density, however, it 
is dominated by stronger divergences. 
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Spatial heterogeneities often give rise to intriguing slow 
dynamics in complex materials, manifested for example by 
broad frequency-dependent relaxation processes in colloidal 
gels which form stress-sustaining networks close to the sol-gel 
transition HI SIS]. Similarly, the presence of differently sized 
proteins, lipids and sugars in the cytoplasm of eukaryotes, 
summarized as cellular crowding, is identified by slow anoma- 
lous transport as its most distinctive fingerprint 
further prominent example are sodium silicates, where the for- 
mation of a space-filling network of channels allows for slow 
diffusion of sodium ions in an arrested host matrix A 
minimal model that encompasses spatial disorder and slow 
dynamics is provided by the Lorentz model [8], i.e., clas- 
sical point particles explore without mutual interaction a d- 
dimensional space in the presence of a frozen array of ran- 
domly distributed (possibly overlapping) hard spherical ob- 
stacles of radius a and concentration n. 

Recently, striking behavior of the velocity auto-correlation 
function (VACF), ip(t) := v~ 2 (v(t) • v(0)), has been re- 
ported for a dense hard-sphere system (d — 3) close to the 
freezing transition jjj]. At intermediate time scales, a regime 
of anti-correlations emerges due to the well-known "rattling" 
of particles in their cages fjjj. The long-time behavior ex- 
hibits an intriguing cross-over scenario from long-living posi- 
tive correlations, ip{t) ~ Af[t~ 3 / 2 , to a high-density regime 
characterized by slowly decaying anti-correlations, ip(t) ~ 
— A' n t~ 5 / 2 . The former corresponds to the celebrated long- 
time anomaly in simple liquids ifTHl . connected to the forma- 
tion of a vortex pattern due to local momentum conservation. 
The mechanism for the latter decay is presumably of totally 
different origin: in an array of immobilized obstacles, the 
dynamics of a tagged particle always remembers its frozen 
cage — a mechanism well-known for the Lorentz model [8]. 

For the Lorentz model itself, there is a long-standing dis- 
crepancy between analytic theory and simulations about the 
manifestation of the long-time tail. The existence of the long- 
time anomaly in the Lorentz model has been predicted within 
a rigorous low-density expansion as ip(t) ~ — A't^ 4 / 2 ^ 1 
for t — > oo J12I Il3f1 . Earlier computer simulations on two- 
dimensional systems HHOE] identified a long-time relax- 
ation of power-law type. At low densities, the expected expo- 



nent was confirmed; the amplitude, however, differed signifi- 
cantly from the theoretical prediction. At intermediate densi- 
ties, the simulation results again suggest power-law behavior, 
which was described phenomenologically by non-universal, 
density-dependent exponents Il5lll7ll . 

The regime of higher obstacle densities poses consider- 
able challenges for theory; sequences of repeated ring colli- 
sions have been accounted for by a self-consistent variational 
repeated-ring theory [18]. Gotze et al. [19] have developed 
a mathematically consistent theory that covers the physics of 
the low-density regime up to the predicted localization transi- 
tion. In particular, they predict a competition between a criti- 
cal power-law relaxation due to the fractal clusters at the per- 
colation transition and the universal long-time tail. A similar 
scenario has been proposed for lattice variants of the Lorentz 
model [20], and there, evidence for a cross-over scenario has 
been reported 112 ill . 

In this Letter, we present high-precision data for the two- 
dimensional overlapping Lorentz model for reduced obstacle 
densities, n* := na d , ranging from the dilute gas, n* = 
0.005, up to the percolation threshold, n* » 0.35907 Q, 
and deep into the localized phase. In particular, we focus on 
the algebraic decay of the VACF at long times which is pre- 
dicted for asymptotically low densities in dimensionless units 
as |UJ] 
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where s — t/r is the mean number of collisions and r _1 — 
2n*v/o- the collision rate 113211 . and v denotes the velocity 
of the particle. In order to observe the universal tails, we 
analyze ip(t) up to times corresponding to several 10 4 col- 
lisions, which implies that a noise level in the correlation 
of the order of 10 _8 (!) is required. Relying on an event- 
oriented algorithm, we have simulated some 10 6 trajectories 
per obstacle density, each trajectory covering 10 6 — 10 8 col- 
lisions. We have calculated ip(t) by directly correlating ve- 
locities and obtain accurate data up to a noise level of 10~ 5 . 
Furthermore, we have measured the mean-square displace- 
ment (MSD), 8r 2 (t) := (AR(i) 2 ), and extracted the diffu- 
sion coefficient D from 5r 2 (t — > 00) ~ 4Dt. An alterna- 
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FIG. 1 : (color online) Velocity autocorrelation function (VACF) for 
the two-dimensional Lorentz model. The minimum is more pro- 
nounced at higher densities. Inset: negative VACF on double loga- 
rithmic scales. Solid lines are fits to the universal long-time tails. The 
universal and the critical power laws are indicated by thick straight 
lines, corresponding to exponents — 2 and 2/z — 2 ~ —1.34. 



tive route to evaluate ip(t) is to perform a numerical second 
derivative of 5r 2 (t) . We have checked that both methods yield 
identical results within statistical errors. The second method 
further suppresses the noise level by several orders of magni- 
tude, up to a factor 10 3 . No smoothening or sophisticated data 
manipulation has been used to obtain the derivatives. 

Results for the VACF are shown in Fig.[T]for the full density 
range in the diffusive phase. On linear scales, the data for the 
lowest density are indistinguishable from the exponential re- 
laxation cxp(— 4t/3r) of the Lorentz-Boltzmann theory. For 
intermediate densities (n* > 0.1), the VACF enters the region 
of anti-correlation already after two collisions. Since the dif- 
fusion coefficient is related to the total area under the VACF 
by a Green-Kubo relation, D = (v 2 /d) J °° \jj{t)At, the ar- 
eas of positive and negative region cancel exactly at densities 
equal and above the percolation threshold, n*. In a double- 
logarithmic representation, the data corroborate power-law 
behavior for time windows covering 1—2 decades or, equiva- 
lently, 2—4 decades in correlations. A gradual increase of the 
density towards n* gives rise to apparent density-dependent 
exponents, at least if correlations below 10~ 4 are ignored. 
Careful inspection of the VACF for n* = 0.20 reveals an in- 
termediate power-law regime as well as an universal long-time 
tail, consistent with the competition of critical and universal 
relaxation predicted by Gotze et al. 01911 . 

As a most sensitive test for the crossover scenario, Fig. [2] 
exhibits the VACF multiplied by the expected power law s 2 
of the universal tail. One infers that s 2 ip(s) saturates at a 
constant in the accessible time window for densities up to 
about 2/3 of the percolation threshold, n* < 0.25. For densi- 
ties n* < 0.1, the constant is approached from above in quali- 
tative agreement with the prediction of Das and Ernst 12311 for 
the sub-leading long-time behavior, 
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FIG. 2: (color online) Rectification of the universal long-time tail 
(a) above and (b) below the percolation threshold. Data points 
are obtained as numerical derivatives of the MSD; open circles for 
n* — 0.20 correspond to directly correlating velocities; extracted 
amplitudes are indicated by solid lines. The dashed line in (b) repre- 
sents the leading correction to the tail, Eq. 



For higher densities, one observes an increase, following 
an apparent density-dependent power law, before the univer- 
sal tail is attained. The time scale where the crossover oc- 
curs shifts to longer times as the percolation threshold is ap- 
proached, confirming the predicted scenario lfl9ll . The density 
n* = 0.10 reaches its asymptotic value remarkably early; this 
is due to a cancellation effect of the sub-leading universal tail 
and the onset of critical slowing down. 

Very close to n*, the VACF exhibits the critical relaxation 
tj){t) ~ t 2 / z ~ 2 , which follows directly from the prediction 
for the MSD, 5r 2 (t) ~ t 2//z . The numerical value of the ex- 
ponent, z w 3.03 B33I1 . coincides with results of simulations 
for diffusion on lattice percolation 02411 . corroborating that the 
transport properties of the Lorentz model in the close vicinity 
of the transition share the same universality class l34Tl . 

Long-time tails originating from power-law distributed exit 
rates of the cul-de-sacs have been predicted even in the lo- 
calized regime |27, 2^1, i.e., for obstacle densities above n*. 
In particular, the VACF should then decay as ip(t) ~ t~ 3 for 
n* > n* and d — 2; a prediction that has not been tested so 
far. Appropriate rectification plots are included in Fig. [2] and 
one infers that the data follow such a power law for one order 
of magnitude in time, i. e., three decades in correlation. 

We have extracted dimensionless amplitudes A of the uni- 



ip(s) ~ -(n/ns 2 ) {1 + 63/16s +...}. 



versal tail as the long-time limit of —s %p(s). In the regime of 
(2) intermediate densities, our data are in semi-quantitative agree- 
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FIG. 3: (color online) Left axis (blue): suppression of the diffusion 
coefficient D/Do with respect to the Boltzmann-Lorentz result Do. 
The dashed line includes the leading low-density correction, the dot- 
ted line corresponds to Eq. l[3}. Right axis (red): reduced amplitude 
A/Ao of the long-time tail from low densities up to the divergence 
at n*. Dashed line: fittoEq. l©, A/A = 1.6n* |n* - n*|~ 91/36 . 



ment with earlier simulations |3 171. see Fig. [3] As has been 
observed before, the values of A are significantly larger than 
the low-density prediction, Ao = n* /it. Since the calculation 
of Ao relies on a perturbative correction to the Boltzmann- 
Lorentz equation 1 13], quantitative agreement requires the dif- 
fusion coefficient D to be sufficiently close to the Boltzmann 
value Do = 3va/16n*. As can be inferred from Fig. [3] this 
criterion is not met even at low densities; a 40% suppression 
of D/Do occurs at n* = 0.10, signalling the onset of sub- 
diffusive motion. The diffusion coefficient for n* < 0.01 is 
in agreement with the first non-analytic correction of the low- 
density expansion lfl2l[l4ll 



Do 
D 



4n* 

— $ = 1- — Inn* - 0.8775 n* + 4.519 (n* Inn*) 2 . (3) 



Although the amplitudes A appear to approach the low- 
density prediction as n* — ► 0, the value for A at n* = 0.005 
still deviates by approximately 25%. 

Close to the percolation threshold, the crossover scenario 
suggests that the amplitude should actually diverge. Match- 
ing the critical relaxation, tp(t) ~ i 2 / z ~ 2 , and the universal 
tail, ip(t) ~ -A(t/r 
t* yields t 2 A ~ t 
crossover of the MSD from anomalous to diffusive transport, 



2/z 



2 , at the divergent crossover time scale 
Assuming that t* also describes the 



2/2 



Dt„, entails the prediction, 

D -2/( Z -2) „ 



T 2 A 



(4) 



as 71* — > n*, where (3 and v are percolation exponents [33]. 
The rapid increase of the amplitudes follows this prediction 
remarkably well; even at n* = 0.1, Eq. (0]i deviates by less 
than 30% from the simulation results, whereas the low-density 
prediction is off by a factor 6, see Fig. [3] 

The VACF or, equivalently, the MSD is only the simplest 
quantity exhibiting anomalous long-time behavior. Devia- 
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FIG. 4: (color online) Top: Critical (solid line) and linear (dashed 
line) divergence of the Burnett coefficient at intermediate densities. 
Bottom: The logarithmic divergence at low densities manifests itself 
as a finite long-time limit of dB(s)/d log(s). 



tions from Fickian diffusion are indicated by a non-vanishing 
(super-)Burnett coefficient, which reads in two dimensions 



Id 

B ^ - 4!d^ 



(AR(t) 2 ) 1 



(5) 



Within a hydrodynamic mode-coupling approach, it has been 
predicted that the Burnett coefficient diverges logarithmically 
in d = 2 E3]. Indeed, d£(i)/dlog(i) saturates at a 
constant for low densities, see Fig. |4] Again close to n*, 
one expects this behavior to be masked by the critical re- 
laxation. Dynamic scaling predicts a power-law divergence 
of the mean-quartic displacement, (AR(t) 4 } ~ t 4 ^ z , where 
z = (2v - + n)/(v - /3/4) « 2.95 < z Q. Then, 
the B{t) increases according to B{t) ~ /j 4 / 2-1 at the critical 
density — consistent with Fig. |U The presence of finite clus- 
ters renders the dynamics spatially heterogeneous, even be- 
low n* . A heterogeneous superposition of Gaussian processes 
yields a linearly divergent Burnett coefficient, B(t) ~ aiD 2 t 
for t — ► 00, with a2 = (4/3)(l/Poo — 1), and Poo denotes the 
fraction of mobile particles. We find semi-quantitative agree- 
ment with this prediction, see Fig. [4] In the dilute limit, the 
prefactor is expected to vanish as (X2D 2 ~ n* . 

Let us briefly comment on the the long-time behavior of the 
VACF for the three-dimensional Lorentz model. Recently, the 
critical properties of the localization transition have been ana- 
lyzed in terms of a scaling Ansatz for the van Hove correlation 
function |25, 2^]. The exponent of the universal tail is larger, 
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FIG. 5: (color online) VACF for the three-dimensional Lorentz 
model; a crossover similar to Fig. Q] from the critical relaxation, 
ip(s) ~ s 2//z_2 , to the universal long-time tail, ip(s) ~ s _5//2 , can 
be observed. 



and the amplitude depends even stronger on the density II 1 3fl 



(37T) 



3/2 



16 



(n*) 2 s~ 5/2 for d = 3 



(6) 



where s = t/r again, and the collision rate is given by t -1 = 
irri*v/a for d = 3. Such a predicted behavior is much more 
difficult to observe, and to the best of our knowledge, only 
rudimentary evidence has been reported for the existence of 
this phenomenon jjlll . Our results (Pig. [5J suggest a similar 
crossover scenario as for d = 2: the universal long-time tail 
with exponent —5/2 is preceded by a critical relaxation with 
exponent 2/z — 2 rj —1.68; the latter covers a growing time 
window upon approaching the localization transition. 

Finally, we emphasize again the universality of the negative 
tail: it relies on the general mechanism that the particle can re- 
turn in its frozen heterogeneous environment by reversing its 
path, thus remembering the presence of free volume |8||. For 
the hard-sphere fluid of Ref. [9], the role of the frozen environ- 
ment is taken by long-lived cages in the high-density regime. 
Hence, the negative tail, tp(t) ~ — A' n t~ 5 / 2 , should emerge in 
an intermediate time window, r -C t <C r a , bounded by the 
time scales for momentum relaxation r and structural relax- 
ation r a ; the slowing down of the latter is also reflected in a 
rapid increase of the viscosity of the fluid, 77 ~ r Q . At even 
larger times, the positive hydrodynamic tail, ip(t) ~ Af\t~ 3 / 2 , 
is expected to follow although its amplitude should vanish 
rapidly, A n ~ r^ 3 / 2 O- 
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